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Abstract: For treatment of the layers below the crust of a neutron star it is use- 
ful to employ a relativistic model involving three independently moving constituents, 
representing superffuid neutrons, superfluid protons, and degenerate negatively charged 
leptons. A Kalb Ramond type formulation is used here to develop such a model for the 
specific purpose of application at the semi macroscopic level characterised by length- 
scales that are long compared with the separation between the highly localised and 
densely packed proton vortices of the Abrikosov type lattice that carries the main part 
of the magnetic flux, but that are short compared with the separation between the 
neutron vortices. 



1 Introduction 

The purpose of this article is to present a concise overview of a class of macroscopic 
relativistic superconducting superfluid models developed |T], Q as a generalisation of pre- 
vious non conducting relativistic superfluid models 0, [|, |5j with a view to applications 
concerning the layers below the crust of a neutron star, which are believed to be well 
described by three constituent superconducting superfluid models of the kind that was 
introduced (as a charged generalisation of the Andreev - Bashkin model||)for a super- 
fluid mixture) by Vardanyan and Sedrakyan^], and that has more recently been further 
developed (though still using a non-relativistic treatment) by Mendell and Lindblom||. 
The three basic ingredients in a description of this kind are, firstly, a condensate of su- 
perfluid neutrons, secondly an independently moving - effectively superconducting - 
condensate of superfluid protons, and thirdly a negatively charged degenerate leptonic 
constituent (consisting mainly of electrons but also including a significant proportion 
of muons) that is of "normal", i.e. non-superfluid, kind. Such a treatment does not 
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include thermal effects (whose inclusion would involve a fourth constituent represent- 
ing entropy) but should nevertheless be applicable as a very good first approximation 
except during a short lived high temperature phase immediately after the birth of the 
neutron star. 

As the relativistic analogue of the kind of phenomenological description introduced 
in a Newtonian context by Bekarevich and Khalatnikov[|9| (as a generalisation of Lan- 
dau's original two-constituent model) it will first be shown how to set up a general cat- 
egory of three-constituent perfectly conducting perfect fluid models of the type that is 
needed, as a preliminary for the more specific developments that follow. This category 
includes, as a specialisation, the case in which the neutronic and the protonic con- 
stituents are both characterised by strictly irrotational behaviour of the kind that is 
relevant in neutron stars on a "mesoscopic" scale, meaning a scale large compared with 
that of the underlying microscopic particle description, but small compared with the 
macroscopic scale separation between the vortex defects within which the superfluid 
comportment is locally violated. Models of this irrotational kind are just a specially 
simple limit within the more general category that is needed for the purpose of treating 
the superconducting superfluid on a "macroscopic" scale, meaning a scale that is large 
compared with the separation between vortices. 

The present discussion will be focussed on an intermediate "semi macroscopic" scale 
meaning a scale that is small compared with the spacing between the superfluid neutron 
vortices (which will be rather widely separated due to the relatively low angular velocity 
of the star, though they contain quite a lot of energy due to their "global" nature) 
but large compared with the Abrikosov lattice spacing between the "local" proton 
vortices, which are expected to be much more numerous in order to carry the rather 
large magnetic fluxes that are thought to be present. 

2 Generic category of 3-constituent superconduct- 
ing superfluid models 

In so far as its contribution to the mass density is concerned, the most important of 
the the three independent constituents under consideration is that of the superfluid 
neutrons, with number current 4- vector n£, say. The second contribution is that of 
the superconducting protons - which make up a small but significant part of the mass 
density - with number current four-vector ng. The third constituent is that of the 
degenerate non-superconducting background of negatively charged leptons - consisting 
mainly of electrons, but including also a certain proportion of muons at the high densi- 
ties under consideration - with a corresponding lepton number number current vector 
say. This negatively charged "normal" (i.e. non superconducting) constituent con- 
tributes only a very small fraction of the mass density, but it nevertheless has a crucially 
important role, not just because the corresponding unit vector u p defined by 



characterises the natural reference frame of rigid corotation in an equilibrium configu- 
ration, but more generally, in so far as electromagnetic effects are concerned, because in 



n£ = n e u p , 
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terms of the electron charge coupling constant e the corresponding total electric current 
4-vector will be given by 

J» = e{n$-n£). (2) 

It will be convenient to express formulae such as this in a condensed notation system 
based on the use of the summation convention for "chemical" indices represented by 
capital Latin letters running over the three relevant values, namely x=n, x=p, x=e. 
Using this convention, the equation (|2]) for the electric current density can be rewritten 
in the consise form 

J p = e*nP, (3) 

where the the charges per neutron, proton, and electron are given respectively by e n = 0, 
e p = e, and e e = — e. 

Since each of the three independent currents involved is conserved, it will be possible 
to use a Kalb-Ramond formulation in which, instead of imposing the three correspond- 
ing conservation laws 

V p < = 0, (4) 

as dynamical equations, they will be obtained as identities by postulating that the 
currents should have the form 

n| = V CT 6f, (5) 

for corresponding antisymmetric gauge bivector fields b pa which are physically defined 
only modulo Kalb Ramond gauge transformations of the form 

K pa ^ K pa + Xd x vpa • (6) 

for arbitrary antisymmetric trivector fields 8 x upcr . 

Since our present treatment will be restricted to the conservative limit in which 
dissipative effects are neglected, the analysis will be assumed to be expressible in terms 
of a variational principle based on a Lagrangian density, in which as usual the repre- 
sentation of the electromagnetic field requires the introduction of a Maxwellian gauge 
1-form A p , in terms of which the gauge invariant electromagnetic field tensor is given 
by 

F pa = 2V [p A a] , (7) 

(using square brackets to indicate index antisymetrisation) . The implementation of 
the Kalb Ramond formulation requires that the set of independent currents n£ be 
supplemented|| by the introduction of a corresponding set of vorticity 2-forms w^ pa 
each of which is characterised both by an algebraic degeneracy condition of the form 

W[ pu w p]a = , (8) 

and by a closure condition of the form 

V[,<] = 0, (9) 

so that each such 2-form w pa is interpretable as a pullback of a prescribed area measure 
on a two-dimensional base space. This means that in terms of suitably chosen local 
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vorticity base space cordinates x 1 ? X 2 the corresponding pair of scalar fields x x , X\ 
induced on the four dimensional spacetime background will specify the vorticities by 
prescriptions of the form w pa = 2x\ p X 2 a } ■ 

In terms of these quantities, a Lagrangian of the appropriate kind will will be ex- 
pressible in the generic form 

£ = A + J p A p + 1 -b° p w* a: (10) 

which consists of a pair of gauge dependent coupling terms preceded by a first term 
A that is a function only of the relevant gauge independent field quantities, which in 
addition, of course to the spacetime metric g pa and the electromagnetic field tensor F pa 
consist of the three independent currents and the three corresponding vorticity 2- 
forms Wp CT . This means that its most general infinitesimal variation will be given by an 
expression of the form 

5A = ft 8n p + iA^ Sw*, + ^H°» 5F pa + ^- 5g pa , (11) 

where the coefficients of the metric variations are not independent of the others but 
must satisfy a Noether type identity of the form 

IT- = IT- = frKsT + + ^-U vp Ff . (12) 

dg P a dg ap 2 2 16vr 

The triplet of covectorial quantities ft is simply interpretable as representing usual 4- 
momenta (per particle) of the neutrons, protons, and leptons. The triplet of rather less 
familiar bivectorial coefficients \° p = — A P<T in this expansion characterises the macro- 
scopic anisotropy arising respectively from the concentration of energy and tension in 
mesoscopic vortices of the neutron and proton superfluids, as a consequence of their 
vorticity quantisation conditions, in the manner discussed in our previous workH on 
the single constituent model. These new four dimensional bivectorial coefficients replace 
the three dimensional (space) vectorial coefficients introduced for a similar purpose in 
a more restricted Newtonian framework by Bekarevich and KhalatnikovQ. Finally the 
bivectorial coefficient 7i pa = —7i pa will be interpretable as an electromagnetic displace- 
ment tensor, in terms of which the total electromagnetic field tensor (0) will be given 
by an expression of the form 

F pa = H pa + 4nM pa , (13) 

in which M. pu is what can be interpreted as the magnetic polarisation tensor. In the 
application that we are considering, this - typically dominant - Abrikosov polarisation 
contribution 47rAi pcr is to be thought of as representing the part of the magnetic field 
confined in the vortices, while the - typically much smaller - remainder 7i pa represents 
the average contribution from the field in between the vortices, which can be expected 
to vanish by the "Meissner effect" in strictly static configurations, but which can be 
expected to acquire a non zero value in rotating configurations due to the London 
mechanism that will be discussed below. 
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In the application of such a generalised Kalb Ramond type variation principle, the 
gauge fields b° p and A p are to be considered as free variables, but n p and are not. 
Each current n p is to be considered as fully determined by the corresponding gauge 
bivector b^ 7 according to the prescription (^) while each vorticity 2 form w* a is to be 
considered as being determined by corresponding freely chosen scalar base coordinate 
pullback fields \\i X 2 ^i which means that the variation of any vorticity 2 form w pa will 
be determined by a corresponding freely chosen displacement vector field £ p according 
to a prescription of the form 5w pv = — 2V[ A1 (uvj p £ p ). 

Subject to these rules, the "diamond" variational integrand 

<>£ = Wgt^Si^gW^C) = 5C + \Cg^ 5g pu . (14) 

needed for the application of the variational principle will be given by 

OC = ( V„ 7T- - \w^)5b^ + /* % + (j p - i- V CT H p °) 5A P + bg pv + V a U% (15) 

in which it is useful to allow for the possibility of varying the background spacetime 
metric g pcr , not only for the purpose of dealing with cases in which one may be concerned 
with General Relativistic gravitational coupling, but even for dealing with cases in 
which one is concerned only with a flat Minkowski background, since, as will be made 
explicit below, the effect of virtual virtuations with respect to the relevant curved 
or flat background can be used for evaluating the relevant "geometric" stress energy 
momentum density tensor T pcr . The coefficients of the current variations are the usual 
gauge dependent total momentum covectors given by 

vr x p = /i x p + e*A p . (16) 

The coefficients of the three independent displacement displacement vector fields 
£ p will be interpretable as the effective force densities acting on the corresponding 
constituent currents. The generic expression for these force densities can be read out 
for the neutrons X=n, protons and X=p and leptons X=e respectively as 

/; = (n; + v,A n >;, (17) 
F P = K + v*vXp > ( 18 ) 

/; = « + v, \rx P ■ (19) 

Although it is of no relevance for the application of the variation principle, it can be 
noted for the record that the current appearing in the final divergence term of (|T5|) will 
be given by 

K ° = tt* 56f - (fof + Af )w*Xl + fabZsr $g,u + ^H™ 5A U . (20) 

An entity of much greater practical interest is the corresponding stress momentum 
energy density tensor, which can be seen to be given by 

T/ = n^l + \?w*, + i«"F w + *tf a , (21) 
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where the generalised pressure function is given by 

v& = A - + b? ( V p K - \w%) . (22) 



The last term in (15) will evidently drop out when we impose the condition of invariance 
with respect to infinitesimal variations of the bivectorial gauge potentials b£ c and b pa 
is imposed, a requirement which can be seen from (|l|) to give field equations of the 
form 

w% = 2V [p 7r x CT] = 2V[^ X CT] + e x F p(T , (23) 

which are evidently equivalent to what in other formulations could be considered just as 
defining relations for the vorticity two-forms. The remaining field equations obtained 
from ( 15f) will consist of 

V a H pu = AttJ p , (24) 
together with the condition that the force density coefficients should all vanish, i.e. 

n = (25) 
for each of the three relevant chemical index values x=n, x=p, x=e. 



3 The semi- macroscopic application. 

To be more specific, we need to specify the scale of application for which the model is 
intended. At a mesoscopic level, meaning on scales large compared with the dimensions 
of individual molecules or Cooper type pairs, but small compared with the intervortex 
spacing, the appropriately specialised model will be of purely fluid type, meaning that 
the function A should not depend on the vorticity forms u> x CT which implies the vanishing 
of the Bekarevich - Khalatnikov coefficients, i.e. the restriction \° p = 0. A further 
restriction that applies at this mesoscopic level is that for the superfluid constituents, 
namely the neutrons and the protons, but not for the degenerate lepton constituent, 
the corresponding vorticities themselves should be zero, i.e. for x^ewe should have 
w* a = 0, which is the integrability condition for the corresponding momenta to have the 
form 27r n p = hV p (f n , 2ir p p = hV p <p p in which the scalars ip n and ip p will be interpretable 
as the phases of underlying bosonic quantum condensates, with periodicity 2n, and 
in which the preceeding factors of 2 have been inserted to allow for the fact that the 
relevant bosons will consist not of single neutrons and protons but of Cooper type pairs 
thereof. 

At a much larger, fully macroscopic scale, involving averaging over large numbers 
of the neutron and proton vortices (that arise as topological defects of the mesoscopic 
phase fields) the neutronic and protonic constituents will be characterised by not just 
by non vanishing effective large scale vorticities w n pa ^ and w p pa ^ 0, but also by non 
vanishing Bekarevich - Khalatnikov coefficients, \° p ^ and \° p ^ 0, so that only the 
degenerate electrons still behave in an effectively fluid manner, in accordance with the 
restriction 

A c CTp = . (26) 
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The purpose of the present article is to focus on an intermediate scale that will 
be refferred to as "semi macroscopic" meaning that it deals with averages over scales 
that are large compared with the spacing between proton vortices, but small compared 
with the spacing between the neutron vortices, which are expected to be relatively 
widely spaced in typical circumstances within neutron stars, whose angular velocities 
are very low as measured by local physical timescales, whereas their magnetic fields are 
typically rather large. On such "semi macroscopic" scales the behaviour of the neutrons 
constituent will not just be of strictly fluid type, meaning that it will be characterised 
by 

K P = , (27) 
but it will also be subject to the mesoscopic superfluidity condition 

27T n p = W p ^ n , (28) 

and hence 

= , (29) 

so that the corresponding dynamical equation, i.e. the requirement (|25|) that the effect 
that the net force density (|T7D should vanish, will be automatically satisfied everywhere 
outside the microscopic cores of the neutron superfluid votices (which are of "global" 
type, meaning that their energy would diverge logarithmically in the absence of the "in- 
fra red" cut off imposed by the presence of neighbouring vortices). However since we 
are considering scales large compared with the separation distance between the much 
more numerous proton vortices (which are of "local" type, meaning that their energy 
density falls off exponentially on a microscopic lengthscale I whose evaluation will be 
discussed below) the proton constituent will be characterised not only by non vanishing 
averaged vorticity, w p pa ^ but also by a non vanishing Bekarevich - Khalatnikov coeffi- 
cient, Ap P 7^ 0. This means that the dynamical requirement that the corresponding net 
force density fll8l) should vanish will provide a rather complicated dynamical equation, 
expressible in the form 

2n; V k //J + en°F ap + w%V v A™ = , (30) 

in which the first term is interpretable as the negative of the Joukovski force density due 
to the "Magnus effect" acting on the proton vortices, the middle term is the Lorentz 
force density representing the effect of the magnetic field on the passing protons, while 
the last term (which was absent in the mesoscopic description) represents the extra force 
density on the fluid due to the effect of the tension of the vortices. As a consequence 
of its "normality" property the leptonic constituent is governed by a dynamical 

equation of the simpler form 

2<V [CT /^+e<F ffp = 0. (31) 

In view of the highly localised nature of the proton vortices, it is reasonable |T0| to 
suppose that their action contribution should be fully determined just by the Abrikosov 
lattice density of such vortices, and that the only independent contribution from the 
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electromagnetic field F pa should be the part provided by the residual - weaker but 
much more widely extended - part of the flux outside the vortex tubes, as given by the 
field 7i pa that is given by the relation (|13|), on the understanding that the polarisation 
contribution <AirM. pa represents the part of the flux attributable to the the vortex tubes. 
This implies that that the gauge independent term A in the action will be decomposable 
in the form 

A = A MV + A F , (32) 

where the macroscopic contribution A MV is required to be functionally independent 
of F^v, so that it is determined just by the vorticity w p and the currents n p , while 
for constency with the variational definition (|ll|) there will be no loss of generality in 
taking the remaining, electromagnetic field dependent, contribution to have the simple 
quadratic form 

a f = -^-n p M ap . (33) 

107T 



which, in the absence of the polarisation contribution AixM. p<J in fll3|) , would reduce just 
to the usual action contribution for an electromagnetic field in vacuum. 

Despite of being much more densely packed than the neutron vortices, the fact that 
(unlike the neutron vortices) the proton vortices are exponentially localised within a 
microscopic confinement radius means that their mutual interactions (unlike those of 
the long range interacting neutron vortices) should remain entirely negligible even for 
extremely high magnetic fields, so that their contribution to the action should be simply 
proportional to their density. This means that it will be possible to make the further 
decomposition 

A M v = A M + A v , (34) 

in which Am is entirely independent of the vorticity w p , so that it depends only on the 
three independent currents n p , while the remainder Ay is just linearly proportional to 
the protonic vorticity density, so that it will be expressible in the form 

A v = X P w p (35) 

where w p is the protonic vorticity magnitude as defined by 



w 



p 



w p P <7 w p pa /2, (36) 



and where, like Am, the coefficient A p depends only on the currents n p . On the basis 
of dimensional considerations - which should be valid provided the Ginzburg landau 
ratio of the London penetration length I that will be discussed below to the relevant 
Pippard correlation length is not too far from the order of unity value that characterises 
the Bogomol'nyi limit [12] - it can be anticipated that, as is confirmed by more detailed 



analysis |13], [14], [H| the Bekarevich - Khalatnikov coefficient A p will have an order of 
magnitude given in terms of the relevant charged particle mass, which in the neutron 
star case under consideration is the proton mass m p (but which in an ordinary metallic 
superconductor would be the electron mass m e ) by by A ~ hn p /m p . 

The appropriate form for the polarisation tensor AnA4 pcr in (|l~3D, can be seen by 
decomposing the vector potential A p as the sum of a gauge dependent contribution 
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proportional to the proton momentum covector 7rP and a gauge independent remainder 



A p in the form 



A p = ^ p + A p , (37) 

from which, by extrior differentiation, one obtains a corresponding decomposition of 
the form 

F pa = ^w p pa + H pcT , (38) 

with 

H pa = 2V[pX] • (39) 
This decomposition has the required form ( |T3"D provided one makes the identification 

M pa = JL< , (40) 
for what will be referred to as the Abrikosov polarisation tensor. 

4 Phenomenological interpretation. 

As discussed in more detail in particular cases 0, |TD[ (correcting earlier workpTS], |T4] in 
which it was overhastily assumed to cancel out) the quantity 4irAi pa defined by (|4~0D 
will be interpretable as representing the part of the magnetic flux confined to the proton 
vortices, whose action contribution will be included in the term Ay given by ([35]), while 
7i pa accounts for the remainder of the flux, which will be distributed over the region 
outside the proton vortices, and whose contribution to the action will be given by fl3~3"|). 
The covector A will be given by 

A p = -- e ^ p , (41) 

so that it will be obtainable from the equation of state function for yP p as derived from 
Am as a linear combination of the form 

A p = A n p + A- + A" , (42) 

in which the terms are proportional respectively to the neutron 4-momentum, the "nor- 
mal" reference state unit vector u p (as specified according to (P by the leptonic current), 
and the (semi macroscopic) electric current J p , so that they will be expressible as 

A% = ~< < , A- = ~ « P , A" = ~^£ 2 J P , (43) 

with proportionality factors that depend (just) on the form of the master function spec- 
fying Am in terms of the relevant currents. In particular the dimensionless parameter a 
would be zero if there were no entrainment, but in view of the effect originally predicted 



by Andreev and Bashkin[y] can be expected[ll] to be of the order of unity, while the ef- 



fective London mass parameter /i L will have a magnitude that is the same as that of the 
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relevant charge carriers, in this case protons, to within a factor comparable with unity, 
from which it differs by an amount that also depends on the entrainment effect. The 
third parameter i is interpretable as the relevant London penetration lengthscale that 
characterises the effective thickness of the individual proton vortices of the Abrikosov 
lattice. If it is assumed, as most authors have done, that the entrainment effect only 
couples the neutrons and protons but does not significantly involve the leptonic back- 
ground (so that the a master function Am of semi separable form0 can be used) then 
it can be estimated that this length scale will be given roughly as a function of the 
effective London mass /z L and the lepton number density n e (which must be very close 
to the proton number density) by 



(44) 



4:Tie 2 n e 

It follows from fl4"3|) that there will be a corresponding decomposition 

n pu = H» a + n- pa + n% , (45) 

with 

W pa = 2V [p Al } , ^ = 2V [p ^l L CT] , H% = 2V [P ^, , (46) 

in which far as the averaged flux is concerned, the main contribution will typically be 
that of the London field meaning the part attributable to the rotation of the 

"normal" (i.e. non-superfluid) negatively charged background. (In an ordinary metallic 
superconductor the analogous London field contribution arises as a well known conse- 
quence of rotation of the positively charged ionic background). In the absence of entrain- 



ment, the neutron vortex contribution 7i. n pcr would vanish. However the expectation |TT 
that the entrainment coefficient a£ will actually be of the order order unity implies that 
although it can be expected to be extremely small outside the immediate neighbour- 
hood of a neutron vortex core (with a confinement radius of the same microscopic order 
of magnitude i as that of a proton vortex) the integrated flux arising from neutron 
vortex contribution TC n pa can be expected to be comparable with that provided by the 
more smoothly distributed (unconfined) London contribution 7i h pa . 

In normal circumstances, the least important term in the sum (|4"5| ) will be what we 
shall refer to as the Meissner residue, meaning the residual contribution Ti.™ arising 
from the semi-macroscopic current J p if any. On the assumption that the lengthscales 
characterising variation of the coefficients fi L and £ are very long compared with the Lon- 
don penetration lengthscale i itself, it can be seen that - except within the microscopic 
defects forming the actual vortex cores where the mesoscopic superfluid description 
breaks down so that w prT is locally non zero - the dominant contribution in the source 
equation (|24| ) will be the part arising from the semi-macroscopic current J M itself, so 
that, as a very good approximation, the source equation will reduce to the well known 
London form 

V CT V CT J"~4 JM > ( 47 ) 

whose homogeneous linear character entails that the only spacially and temporally 
uniform solution is that for which the current J p simply vanishes. What this implies 
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is that after any initial high frequency osillations that may have been present have had 
time to radiate away or otherwise be dissipated, the medium will tend to settle towards 
a state in which the current actually is zero outside a very small radius of order I 
surrounding each individual vortex, so that its average ( J p ) over scales large compared 
with i will also tend to vanish, 

(J')~0. (48) 

The same conclusion therefore applies to the corresponding residual Meissner field con- 
tribution, W™, which will end up in a state such that 

(H%) * . (49) 

This last result is interpretable as a generalisation to "type II" (London - Abrikosov) 
superconductors of the Meissner effect that was originally observed in laboratory exam- 
ples of "type I" , meaning cases in which the Ginzburg Landau ratio of the penetration 
lengthscale i to the relevant Pippard correlation lenghthscale is so small that, instead 
of condensing into an Abrikosov vortex lattice, the magnetic flux tends to be entirely 
expelled into domains where the superconductivity breaks down. In a type I situation 
the polarisation M. pa associated with the Abrikosov lattice will be absent, so there will 
be no distinction between the total flux F pa and the contribution 7i p(J as defined here, 
which means that in this experimentally more familiar case, expulsion of 7i pv is equiv- 
alent to complete expulsion of F pcr . On the other hand in the type II case, although 
there will be the same tendency to expulsion of 7i pcr , this will not entail the complete 
expulsion of F pa because the Abrikosov polarisation contribution M. pv will still remain. 

Both in the type II and - as originally remarked by London - in the type I case, 
the tendency for 7i pv to vanish will be partially thwarted in a rotating background, for 
which the London contribution TC^ will still remain even after the residual Meissner 
contribution TC^ has been dissipated in accordance with (f49|). In the usual laboratory 
applications, whether of type I or type II, this London contribution is all that will 
remain, but in the neutron star case there will also be the neutron vortex contribution 
TOL . If the gradient of the (weakly density dependant) coefficient a£ is not entirely 
negligible, then as well as having a dominant part that is of magnetic character, may 
also include a small contribution to the electric displacement vector D p as defined with 
respect to the "normal" leptonic background frame by the decomposition 

n pa = H pa + 2u [p D a] , D p = H pa u° . (50) 

Since the alignement of the neutron momentum covector will not on average be very 
different from that of the background frame, it can be seen that what is to be expected 
is that under typical equilibrium conditions the macroscopically averaged value of the 
first contribution in ( [45]) will be given by 

(H%) ~ 2u [p (D^ + (H^) , (51) 

in which the main contribution is the purely magnetic part given by 

(52) 
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while the - typically very small - electric part will be given by 

(fl^-^AX), (53) 

where the relevant neutron Fermi energy parameter is given by /x n = —u a fi^. 

The analogous macroscopic average of the London contribution can be instructively 
formulated in terms of the "normal" background's acceleration tensor ii p and rotation 
tensor Q pcr (whose magnitude Q = ^Vtp a Vt pa /2 is the local angular velocity) as defined 
by 

V[ p u a ] = Q pa - u [p u a] , u p = u a V a u p . (54) 

The ensuing result (correcting a misplaced factor of 2 in the preceeding presentation^]) 
is expressible as 

(H^) ~ 2u lp {Dl ] ) + (H- a ) , (55) 
in which the main contribution is the magnetic part which can be seen to be given by 

(H^) = -j>HV> (56) 

As originally observed by London, this is proportional to the angular velocity Q, with 
a proportionality factor /i L that in the neutron star application will be given roughly 
(but due to the "entrainment" effect not exactly) by the proton mass (whereas in the 
ordinary metallic superconductors originally envisaged by London it is given by the 
electron mass). As well as this well known magnetic contribution, there will also be 
a corresponding, but typically much less important, electric displacement contribution 
given by 

(I>) ~ -(/^ + V^ L >. (57) 

H e 

(This electric displacement field is usually ignored in discussions of laboratory applica- 
tions, but even in the ideally simplified case of a motionless incompressible sample that 
is postulated to be strictly homogeneous so that the gradient term would be absent, 
a small residual electric field of this type would still be needed to balance the effect 
on the conducting particles - which in that case would be electrons - of the ordinary 
terrestrial gravitation field.) 

I wish to thank Silvano Bonazzola, David Langlois, and Reinhard Prix for many 
helpful discussions, and I particularly wish to express my appreciation to Isaac Kha- 
latnikov for instructive conversations on numerous occasions since my interest in this 
subject was originally inspired by his classic textbook [[T5[]. 
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